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ABSTRACT 


Lidar  (laser  radar)  data  obtained  at  Hamilton  AFB,  California, 
under  conditions  of  low  ceiling  and  visibility  are  analyzed  by  hand  and 
by  electronic  computer  to  explore  the  operational  utility  of  lidar  in 
determining  cloud  ceiling  and  visibility  for  aircraft  landing  operations. 
Hand  analyses  of  the  data  show  the  ability  of  the  lidar  to  describe 
the  spatial  configuration  of  the  low-cloud  structure  along  the  landing- 
approach  path.  The  problems  inherent  in  evaluating  lidar  observations 
are  discussed,  and  initial  approaches  to  quantitative  solutions  by 
computer  are  presented.  It  is  demonstrated  that  operationally  useful 
information  on  the  ceiling  and  visibility  conditions  contained  in  the 
hand  analyses  can  be  presented  by  digitizing  the  lidar  data  and  sub¬ 
jecting  these  data  to  computer  analysis. 
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I  I  NTRODl'CT  I  ON 


On  8  and  9  January  1968,  the  SRI  Mark  V  pulsed  ruby  lidar  was 
activated  at  Hamilton  AFB,  California,  in  order  to  explore  the  opera¬ 
tional  utility  of  the  lidar  in  cloud  ceiling  and  visibility  determina¬ 
tion  for  aircraft  landing  operations.  Because  of  a  unique  location  of 
the  airfield  on  the  western  edge  of  San  Pablo  Bay  (see  Fig.  1),  the 
aircraft  landing  operations  are  confronted  with  a  difficult  meteoro¬ 
logical  problem.  San  Pablo  Bay  with  its  relatively  cold  water  is  a 
notorious  source  of  fog  and  low  stratus,  particularly  in  winter.  The 
landing-approach  glide  path  begins  over  the  marshes  and  open  water  of 
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FIG.  1  LOCATION  OF  SRI  LIDAR  AT  HAMILTON  AFB  RELATIVE 
TO  TOUCHDOWN  POINT  AND  TO  SAN  PABLO  BAY 
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the  Bay,  an  area  cl.-racterizea  by  low  ceiling  and  visibility,  and  extends 
six  miles,  with  a  2.5-degree  slope,  to  the  point  of  touchdown.  The  Air 
Weather  Service  observing  station  and  a  rotat ing-beam  ceilometer  are 
located  near  touchdown — i.e.,  as  near  as  possible  to  the  area  of  deterio¬ 
rating  ceiling  and  visibility  conditions  (see  Fig.  1).  However,  even 
at  this  forward  location,  ceiling  and  visibility  conditions  are  often 
quite  different  from  those  encountered  by  incoming  aircraft  over  the 
marshes  and  open  water. 

At  the  request  of  the  base  commander  and  the  weather  squadron 
commander,  SRI  considered  the  application  of  lidar  to  their  special 
problem  of  remotely  measuring  ceiling  and  visibility  along  the  landing 
approach  path.  The  field  experiment  could  not  involve  a  large-scale 
effort,  and  therefore  the  collected  data  are  incomplete  in  many  respects. 
Consequently,  analyses  and  discussions  of  the  data  presented  in  this 
report  must  be  considered  exploratory. 

The  manner  in  which  the  lidar  data  were  collected  and  analyzed  is 
discussed  in  Sec.  II.  Results  of  the  hand  analyses  are  presented  in 
Sec.  III.  Computer  techniques  to  numerically  process  lidar  data  related 
to  cloud  ceiling  and  visibility  are  discussed  in  Sec.  IV. 

The  potential  application  of  lidar  to  remote  measurements  of 
ceiling  conditions  is  considered  promising,  particularly  in  view  of 
the  future  availability  of  higher-performance  lidars. 
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II  METHOD  OF  OBSERVATION  AND  DATA  ANALYSIS 


The  SRI  Mark  V  pulsed  ruby  lidar,  details  of  which  are  given  in 
Table  I,  was  transported  to  Hamilton  AFB  and  was  placed  next  to  the 


Table  I 

CHARACTERISTICS  OF  SRI  MARK  V  RUBY  LIDAR 


Transmitter 

Laser 

6X  3/8-inch  ruby  crystal,  Brewster-angle 
one  end,  planar  one  end,  uncoated 

Q-switch 

Rotating  prism  and  saturable  dye 

Wavelength 

6943  A 

Pulse  length 

ia  ns 

Peak  power  output 

18  MW 

Pulse  energy  output 

0.27  Joule 

npt ics 

6-inch  Newtonian  reflector  telescope 

Beamwidth 

Approximately  0.3  mrad 

PRF 

Two  per  minute 

Receiver 

Detector  photomultiplier 

14-stage  RCA  type  7265 

Optics 

6- inch  Newtonian  reflector  telescope 
with  adjustable  field  stop 

Field  of  view 

0.2  to  0.9  mrad 

Bandpass 

Approximately  17  A 
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*iMthor  observing  station  and  the  rota t mg -beam  ceilometer.  Data 
related  to  cloud  ceiling  and  visibility  were  collected  by  firing  at  a 
pulse  rate  frequency  ol  one  or  two  pulses  per  minute  out  across  the 
Bay.  parallel  to  the  aircraft  glide  path  (see  Fig.  1).  The  elevation 
angle  of  the  direction  of  firing  was  varied  from  zero  to  65  degrees 
without  change  of  azimuth.  All  data  were  obtained  under  the  actual 
weather  conditions  that  create  the  operational  problems.  The  weather 
conditions  that  prevailed  during  the  period  of  lidar  operation  are 
given  in  lable  II.  Lidar  observations  were  made  during  the  afternoon 
and  evening  ol  8  January  and  during  the  forenoon  of  9  January. 


Table  II 

WEATHER  CONDITIONS  DURING  LIDAR  OBSERVATIONS  AT  HAMILTON  AFB,  CALIFORNIA 


Weather  Conditions 

8  January  1968 

9  January  1968 

Cloud  ceiling 

700-800  ft 

■100-500  ft 

Prevailing  visibility 

1-1  2  to  2-12  mi 

1  1  to  1/2  mi 

Obstruction  to  visibility 

Fog,  occasional 
light  rain 

Fog,  light  drizzle 

Temperature 

37  -38F 

39  -10  F 

Each  lidai  firing  was  separately  recorded  by  photographing  the 
trace  of  received  signal  power  vs.  slant  range  (often  called  the  back- 
scatter  signature)  as  it  appeared  on  the  oscilloscope  (A  scope'). 
Interpretation  of  the  recorded  data  is  best  explained  with  the  aid  of 
the  lidar  equation  which  can  be  written  in  the  form 

CT  Ar  r  rR  ^ 

Pr  ([i)  =  Pt  T~  S180(R)T-(R)  cxp  "  2  '  CT(R')dR'j  (1) 

K“  ~  L  v  o  J 


P  -  Power  collected  (at  a  given  instant)  by  the  primary 

receiver  optics  from  atmospheric  backscatter  of  laser 
energy 
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R  =  Range 


=  Power  transmitted  into  the  atmosphere 

c  =  Velocity  of  light 

t  =  Pulse  duration  (seconds) 

A  =  Effective  receiver  area 
r 

B  (li)  -  Volume  backseat  ter  coefficient 
1  oO 

T  (R)  Beam  convergence  factor 

3(10  =  Volume  extinction  coefficient 

The  above  lormulation  assumes  a  constant  energy  density  across 

the  beam,  randomly  distributed  scatterers  within  the  effective  scattering 

volume,  and  Bouquer's  law  of  attenuation.  When  the  beam  convergence 

factor,  T  .  is  excluded,  Eq .  (1)  is  valid  only  at  ranges  beyond  the 

point  where  the  diverging  transmitted  beam  of  the  non-coaxial  lidar 

system  is  fully  encompassed  by  the  diverging  receiver  field  of  view. 

By  including  T  ,  Eq.  (1)  can  describe  the  behavior  of  the  lidar  return 
c 

signal  at  close-in  range.  varies  between  0  (before  beam  interception) 

and  1  (at  full  beam  convergence)  and  is  a  function  of  the  specific  lidar 
system  used. 

Figure  2  shows  examples  of  lidar  data  obtained  at  three  elevation 

angles.  The  received  signal  power  (related  to  P  )  on  a  logarithmic 

r 

scale  is  recorded  vs.  slant  range  (R)  on  a  linear  scale.  The  use  of 
a  logarithmic  video  amplifier  in  the  receiver  is  almost  essential  in 
order  to  compress  the  wide  dynamic  range  of  the  detector  (typically 
four  or  more  orders  of  magnitude)  and  to  enable  the  received  signal  to 
be  displayed  on  a  single  oscilloscope  trace  without  loss  of  detailed 
information.  The  sweep  speed  of  the  oscilloscope  in  Fig.  2(b)  and 
2(c)  is  twice  as  fast  (1  microsecond  per  division  with  a  maximum  re¬ 
corded  range  of  1.5  km)  as  that  in  Fig.  2(a)  (2  microseconds  per  division 
with  a  maximum  recorded  range  of  3.0  km)  in  order  to  show  the  significant 
lidar  echoes  detected  at  the  higher  elevation  angles  in  more  detail. 

At  close-in  ranges  (up  to  about  100  m)  where  the  beam  convergence  factor 
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complete  "scan"  from  the  horizontal  to  the  near  vertical  are  analyzed 
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and  related  to  the  low  cloud  structure  and  the  cloud  ceiling  as  measured 
by  the  ceilometer. 

Recorded  data  of  received  signal  power  vs.  range  obtained  along 
the  horizontal  line  of  sight  can  be  related  to  horizontal  visibility 
while  those  obtained  at  low  elevation  angles  [Fig.  2(a)]  can  be  related 
to  slant  visibility.  However,  for  correct  application  to  visibility, 
the  lidar  data  need  to  be  processed  numerically  to  extract  the  volume 
extinction  coefficient,  r(R),  which  is  fundamental  to  the  determination 
of  visibility.  Because  of  the  limited  scope  of  the  Hamilton  AFB  experi¬ 
ment,  the  transfer  characteristics  of  the  receiver  optics  (including 
photomultiplier,  logarithmic  amplifier,  and  oscilloscope)  could  not  be 
determined  with  the  required  accuracy.  Therefore,  values  of  received 
signal  power  could  not  be  reliably  transformed  into  absolute  values  of 

P  (R).  Furthermore,  examination  of  the  lidar  data  after  the  experiment 
r 

revealed  that  the  very  large  backscatter  signal  from  the  fog  at  close- 
in  range  tended  to  saturate  the  receiver.  This  receiver  saturation, 
in  turn,  caused  a  time-dependent  discharge  of  photons  ("after  pulsing") 
which  invalidates  any  computations  based  on  the  slope  of  the  backscatter 
curves  obtained  at  low  elevation.  In  view  of  these  difficulties  no 
worthwhile  analysis  of  lidar  data  in  terms  of  horizontal  and/or  slant 
visibility  is  presented. 
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Ill  LIDAR  OBSERVATIONS  RELATED  TO  CUR'D  CEILING 

A.  Hand  Analysis  of  Available  Data 

Figures  3,  4,  and  5  present  the  spatial  distribution  of  the  low- 
level  cloud  layers  as  observed  with  the  lidar  in  the  direction  of  San 
Pablo  Bay  during  three  separate  time  periods  on  8  January  1968.  The 
spatial  distributions  are  obtained  by  firing  the  lidar  at  successive 
elevation  angles  from  a  maximum  of  65  degrees  down  to  the  horizontal 
and  analyzing  along  each  direction  of  firing  all  echoes  related  to 
cloud  layers  and  to  levels  of  optical  density  change.  Layers  are 
indicated  by  solid  bars  and  levels  by  solid  triangles.  Representative 
samples  of  individual  lidar  shots  accompany  each  figure.  Also  indicated 
in  each  I igure  is  the  available  surface-weather  observation  for  the 
time  closest  to  the  period  of  the  lidar  data.  Observed  cloud-ceiling 
height  (in  hundreds  of  feet)  is  given  by  the  ceilometer.  Observed 
horizontal  visibility  is  the  so-called  prevailing  visibility  defined  as 
’’the  greatest  visibility  that  is  attained  or  surpassed  throughout  at 
least  half  ol  the  horizon  circle,  not  necessarily  continuous.” 

Figures  3  through  5  show  that  under  the  prevailing  weather  conditions 
the  lidar  is  capable  of  describing  the  low-cloud  structure  from  the  obser¬ 
vation  site  out  to  a  distance  of  0.8  to  1.2  km  (1/2  to  3/4  mi).  This 
relatively  limited  range  is  due  to  the  large  attenuation  of  the  lidar 
pulse  energy  along  the  slant  paths  through  the  fog  at  low  elevation. 

During  the  evening,  when  the  reported  horizontal  visibility  increases 
from  1-1  2  to  2-12  miles,  the  range  of  lidar  cloud-detection  reaches 
the  maximum  of  3  4  mile.  The  cloud-ceiling  height  measured  at  the 
observation  site  by  the  rota t ing- beam  ceilometer  corresponds  to  the 
height  of  the  lowest  layer  of  lidar  echoes  obtained  at  high  elevation 
angle . 

Figure  3  shows  a  marked  difference  in  the  cloud  structure  between 
elevation  angles  larger  and  less  than  about  35  degrees.  This  difference 
probably  arises  from  the  fact  that  at  high  elevation  angle  the  low-cloud 
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echoes  are  detected  at  close  range,  with  a  subsequent  loss  of  resolution 
in  the  recorded  lidar  trace  (see  sample  of  individual  lidar  shot  at 
65  degrees  elevation).  Under  such  conditions  a  change  in  the  sweep 
speed  of  the  "A  scope''  must  be  used  to  "magnify'  the  lidar  trace. 

Figures  4  and  5  show  that  higher-level  cloud  layers  are  detected 
only  at  the  high  elevation  angles,  where  the  path  length  through  the 
lower  clouds  and  the*  fog  is  minimum.  The  ceilometer  detected  higher 
cloud  layers  only  during  the  time  period  of  Fig.  5  when  the  ceiling 
was  broken. 

Figures  6  and  7  illustrate  the  capability  of  the  SRI  ruby  lidar 
to  remotely  monitor  cloud-ceiling  variations  at  a  point  distant  from 


FIG.  6  TIME  VARIATION  OF  LOW-CLOUD  STRUCTURE  MONITORED  BY  LIDAR  AT 
1  2  TO  3  4  MILE  FROM  THE  LIDAR  SITE  AT  HAMILTON  AFB  DURING 
AFTERNOON  OF  8  JANUARY  1968 


the  observation  site.  Figure  6  shows  a  12-minute  time  section  of  the 
low- level  cloud  structure  over  the  Bay  at  points  approximately  1/2  to 
3  1  mile  from  the  lidar  site  (1!  degrees  elevation  angle)  during  the 
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afternoon  of  8  January.  The  time  period  corresponds  to  that  of  Fig.  3. 

The  levels  of  signal  Increase,  peak  signal  return,  and  signal  decrease 
monitored  at  time  intervals  of  about  30  seconds  are  indicated  and  joined 
by  straight-line  segments  to  portray  the  time  variation  of  a  100-meter- 
thick  low-cloud  layer.  Practically  no  change  with  time  is  evident. 
Conditions  at  the  remote  locations  are  nearly  identical  to  those  "over¬ 
head"  at  the  observation  site,  the  height  of  the  peak  signal  return 
being  nearly  identical  to  t  lie  ceiling  height  given  by  the  ceilometer. 
Figure  7  shows  the  low-cloud  configuration  as  recorded  with  the  lidar 
at  a  point  about  1  -1  mile  I  rom  the  observation  site  (30  degrees  elevation 
angle)  during  a  time  period  of  3  to  1  minutes.  The  data  are  part  of  the 
series  analyzed  in  Fig.  5  and  were  obtained  on  the  evening  of  8  January. 
The  surface  weather  observation  before  and  after  the  period  of  lidar 
observations  is  indicated.  The  lowest  cloud  layer  detected  by  the  lidar 
(indicated  in  Fig.  7  by  arrows)  closely  corresponds  to  the  800-ft  (244  m) 
ceiling  measured  by  the  ceilometer.  On  two  occasions  (22:33.5  LST  and 
22:31  LST)  the  lidar  data  show  the  upper  (1500  ft)  cloud  layer  detected 
by  the  ceilometer  at  21:58  LST.  Thus,  the  density  variations  in  the 
lower  elouds  apparent  from  the  lidar  data  are  recorded  by  the  ceilometer 
a  1  sci. 

During  the  morning  of  9  January,  very  low  ceiling  and  visibility 

prevailed.  Up  to  11:57  LST  the  sky  cover  remained  overcast,  with  the 

cloud  base  varying  between  400  ft  (122  m)  and  500  ft  (152  m).  The 

prevailing  horizontal  visibility  remained  around  12  mile  (0.8  km)  in 

light  drizzle  and  fog  (see  Table  II).  At  11:57  LST,  clouds  became 

broken.  Figure  8  illustrates  some  typical  lidar  backscatter  profiles 

obtained  throughout  the  morning  at  30  degrees  elevation.  The  rapid 

decrease  in  return  signal  amplitude  observed  at  a  height  near  150  m 

(192  ft)  at  11:33  LST  must  be  interpreted  as  due  to  a  rapid  increase 

in  optical  density,  since  this  height  corresponds  to  that  of  the  cloud 

ceiling  measured  by  the  ceilometer  at  the  lidar  site.  Furthermore, 

* 

at  11:50  LST  an  airborne  weather  observer  reported  a  400-ft  (140-m) 

Capt .  R.  H.  Hcdenberg,  Commander,  Det.  9,  35th  Weather  SQ. ,  Hamilton  AFB. 
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FIG.  8  SAMPLES  OF  LIDAR  DATA  OBTAINED  AT  30  DEGREES  ELEVATION 
AT  HAMILTON  AFB  ON  9  JANUARY  1968 


cloud  ceiling  at  1  mile  from  touchdown  and  indicated  that  clouds  ex¬ 
tended  vertically  to  at  least  2000  ft  (610  m).  Near  12  o'clock  noon, 
however,  this  dense  lower  layer  became  transparent,  as  evident  by  the 
appearance  of  a  cloud  echo  near  360  m  in  the  lidar  profile  of  11:58  LST. 
The  higher- level  cloud  layer  was  also  detected  by  the  ceilometer  at 
the  lidar  site  and  was  firs'  reported  in  the  surface  weather  observation 
at  11:57  LST. 

Figure  9  shows  the  spatial  extent  of  the  significant  layers  and 
levels  detected  by  the  lidar  shortly  after  noon.  The  lowest-alt itude 
signals  correspond  to  the  400-ft  broken-cloud  layer  indicated  by  the 
ceilometer.  Agreement  between  ceilometer  and  lidar  data  extends  to 
1200  ft. 
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12:11  LST  (Solid  Bars:  Cloud  Layers;  Solid  Triangles:  Levels  of  Optical  Density  Change) 


B. 


Discussion 


The  data  samples  analyzed  in  Figs.  3  through  9  demonstrate  the 
advantages  and  limitations  of  the  lidar  in  the  description  of  the  low- 
level  cloud  structure  under  adverse  weather  conditions.  Because  of 
its  pulsed  nature,  the  lidar  can  provide  a  vertical  density  profile 
through  the  low-level  clouds  and,  furthermore,  because  of  its  ability 
to  operate  under  variable  elevation  and  azimuth,  can  accurately  describe 
the  spatial  distribution  of  the  cloud  ceiling.  A  most  useful  feature  of 
the  lidar  is  its  ability  to  monitor  the  cloud-ceiling  conditions  at  a 
location  that  is  remotely  situated  with  respect  to  the  lidar  site. 

Under  the  weather  conditions  that  prevailed  during  the  Hamilton  AFB 
experiment,  the  lidar  could  accurately  monitor  the  cloud  ceiling  at  a 
point  1,2  to  3  4  mile  from  touchdown  point. 

A  replacement  of  the  ceilometer  by  a  lidar  would  greatly  expand 
the  range  of  information  on  cloud-ceiling  height  and  low-cloud  structure 
that  could  be  obtained.  On  the  other  hand,  the  lidar  data  show  a  high 
spatial  correlation  with  the  spot  information  provided  by  the  ceilometer. 
Under  the  conditions  of  fog  and  drizzle  that  prevailed  during  the  ex¬ 
periment,  the  ceilometer  reading  appeared  equally  valid  for  locations 
up  to  1/2  to  3  4  mile  away.  With  the  experimental  lidar  used,  it  was 
disappointing  that  the  maximum  range  at  which  the  low  clouds  could  be 
accurately  described  remained  well  below  the  1-mile  mark.  During  the 
conditions  of  fog  and  light  rain  or  light  drizzle  that  prevailed  during 
the  experiment,  the  long  slant  paths  at  elevation  angles  below  10  to  15 
degrees  rapidly  attenuated  the  ruby  lidar  pulse  energy  down  to  the  noise 
level.  Improvement  in  the  range  of  cloud-mapping  capability  can  be 
expected  in  the  future  by  further  narrowing  the  transmitted  lidar  beam 
"•nd  by  introducing  higher-power  output. 
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IV  COMPUTER  TECHNIQUES  FOR  ANALYZING  LI  DAR  DATA 
RELATED  TO  CLOUD  CEILING  AND  VISIBILITY* 

This  section  describes  methods  of  producing  computer-analyzed 
cross  sections  of  atmospheric  optical  parameters  from  an  input  of 
digitized  lidar  data.  Because  of  the  above-mentioned  instrumental 
uncertainties  with  respect  to  receiver  transfer  characteristics  and 
overloading,  results  of  the  numerical  analyses  must  be  evaluated  in 
terms  of  relative  variations  rather  than  absolute  values — i.e.,  appli¬ 
cation  is  to  cloud  ceiling  rather  than  to  visibility. 

A •  Computer  Analysis  of  Data  Sample 

In  order  to  digitize  the  recorded  lidar  data,  Polaroid  prints  of 
the  oscilloscope-displayed  lidar  signatures,  such  as  shown  in  Fig.  2, 
are  projected  with  a  magnification  of  5  power  onto  a  screen  of  a  scaling 
machine  (digitizer)  that  records  on  IBM  cards  the  (X,  Y)  coordinates  of 
a  movable  crosshair.  The  operator  records  (X,  Y)  values  of  all  in¬ 
flection  points  on  a  lidar  recorded  signature  that,  when  joined  by 
straight  lines,  provide  a  digitized  representation  of  the  signature. 

The  (X,  Y)  representation  of  each  lidar  signature  thus  obtained  consti¬ 
tutes  the  input  for  a  computer  program,  which,  by  use  of  the  digitizer 
calibration  and  the  beam  elevation  angle,  converts  the  (X,  Y)  values 
into  a  matrix  of  oscilloscope  deflection  (output  voltage)  of  the  data 
points  and  a  matrix  of  position  of  the  data  points.  Such  output  data 
from  multiple  laser  firings  obtained  as  the  lidar  scans  from  the  hori¬ 
zontal  to  the  near  vertical  are  then  fed  into  a  grid-point-analysis 
program.  This  program  assigns  a  value  of  the  dependent  variable  (in  this 


*Pa j i  ol  the  work  reported  on  in  this  section  was  supported  by  Stanford 
Research  Institute  (see  Acknowledgments)  and  SRI  Project  7165.  (See 
Ref.  1.  References  are  listed  at  the  end  of  the  report). 

* 

Developed  by  Mr.  R.  L.  Mancuso  of  the  ncrophysics  Laboratory. 
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case,  oscilloscope  deflection  in  volts)  to  each  grid  point  of  a  pre¬ 
determined  grid  in  the  vertical  plane  of  observation.  Values  to  in¬ 
dividual  grid  points  are  assigned  on  the  basis  of  the  five  data  points 
nearest  to  the  grid  point  with  the  aid  of  a  distance-weighting  factor 
as  well  as  a  vector  weight  defined  in  the  radial  direction  (which  gives 
more  weigtit  to  data  from  a  single  lidar  profile)  or  in  the  horizontal 
direction  (which  gives  weight  to  existing  conditions  of  horizontal 
st  ra t i f icat ion ) . 


Figure  10(a)  shows  the  results  of  this  type  of  computer  analysis 
using  the  lidar  data  observed  during  the  evening  of  8  January  1968 
(22:22-22:31  LST ) .  The  analysis  is  prepared  from  the  data  of  18  in¬ 
dividual  lidar  profiles  obtained  at  elevation  angles  ranging  from  zero 
to  60  degrees.  In  Fig.  10  and  in  all  subsequent  figures,  coordinate 
axes  and  isopleths  are  drawn  by  hand.  In  order  to  supply  as  much  input 
data  as  possible  to  the  grid-point-analysis  program,  data  points  inter¬ 
mediate  to  the  inflections  were  obtained  by  linear  interpolation  within 
the  computer  program.  Minimum  distance  between  signature  data  points 
was  25  m.  The  oscilloscope-deflection  analysis  of  Fig.  10(a)  provides 
a  field  that  is  representative  of  atmospheric  scattering  activity.  The 
layer  of  maximum  scattering  activity  at  a  height  of  250  to  275  m  is 
related  to  the  level  of  low  clouds  shown  in  the  hand  analysis  of  the 
same  date  in  Fig.  5. 


It  is  possible  to  transform  the  voltage  f.eld  of  Fig.  10(a)  into 
a  parameter  field  more  representative  of  atmospheric  conditions  and 
perhaps,  in  the  case  of  ideal  data,  into  a  field  of  extinction  coef¬ 
ficients  from  which  the  meteorological  range  or  the  "visibility  along 
any  given  path  (such  as  along  an  aircraft  landing  approach  path)  may¬ 
be  inferred.  Referring  to  Eq .  (1),  the  lidar  equation,  a  convenient 
range-corrected  quantity,  in  dB  notation,  can  be  defined  as: 


5(H)  =  10  log 


P(R)R 


P(R  )R 
o  o 


10  log 


B180(R)Ta(R)Tc(R) 

6180(Ro)Ta2(Ro)Tc(Rc,) 


(2) 
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FIG.  10  COMPUTER  PRINTOUT  FOR  GRID-POINT  ANALYSIS  OF  (o)  OSCILLOSCOPE-DEFLECTION  IN  RELATIVE 
VOLTAGE  AND  (b)  S-FUNCTION  IN  RELATIVE  dB,  IN  THE  VERTICAL  PLANE  OF  LlDAR  OBSERVATION 
(HAMILTON  AFB,  8  JANUARY,  1968,  22:22  -  22:31  LST) 


whore  Rq  is  a  relcrence  range  on  each  lidar  trace  and 


R 

r  ,  , 

T  (R)  -  exp  -  |  d(R  )  dR 
v  o 


is  the  one-way  atmospheric  transmission  over  the  range  R.  All  other 
terms  are  as  defined  previously.  Again  using  the  lidar  data  observed 
during  the  evening  ol  8  January,  a  grid-point  analysis  of  S(R)  was 
obtained  in  the  following  manner.  From  a  relative  calibration  of  the 
nonlinear  response  of  the  Mark  V  ruby  lidar  receiver  (including  photo¬ 
multiplier,  logarithmic  amplifier  and  oscilloscope)  estimated  to  be 
valid  for  the  conditions  of  the  Hamilton  AFB  experiment,  a  polynomial 
operator  was  defined  by  means  of  which  the  computer  transforms  the 
output-vol tage  field  of  Fig.  10(a)  into  a  field  of  relative  light 
IP  (R)J  incident  on  the  primary  receiver  optics.  Using  the  first  data 

point  on  each  ol  the  18  lidar  observations  as  the  reference  range  R  , 

o 

Fig.  10(b)  shows  the  computer  output  of  the  S-funetion  field  thus 
obtained.  Large  increases  and  decreases  with  height  correspond  to  cloud 
layers.  A  comparison  with  the  qualitative  analysis  of  the  same  data 
shown  in  fig.  3  clearly  indicates  the  practicability  of  this  computer- 
produced  parameter  field  for  cloud-ceiling  determination. 

The  derivative  of  8(K)  over  the  region  of  full  beam  convergence 
(T  =  1)  .'lay  be  expressed  as: 


dS 

dR 


4.34 


1  dP 

I  dR 


8.7c 


(3) 


where  the  subscript  180  and  the  range-dependence  notation  (R)  have  been 
omitted.  A  closed-form  solution  for  the  optical  parameters  (3  and  C 
requires  additional  information  on  the  scattering  properties  of  the 
atmosphere  or  on  a  backsca t t er-ext inct ion  relation.  When  the  backscatter 
coefficient,  p.  is  independent  of  range,  as  is  the  case  in  a  homogeneous 
scattering  medium,  the  extinction  coefficient  between  deflection  points 
of  8(R)  is  given  by 


1  LS 
8,7  iR 


(4) 
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Figure  11(a)  presents  a  computer  output  of  the  grid-point  analysis 
ot  c  over  the  complete  field  of  Fig.  10.  Negative  values  (shaded)  and 
large  positive  values  (unshaded)  represent  areas  where  the  received 
lidar  signal  sharply  increases  and  decreases,  respectively.  Beyond 
the  range  of  lull  beam  convergence,  these  areas  correspond  to  areas  of 
non-constant  backseat  ter  and  allow  inferences  as  to  the  levels  of  low- 
cloud  layers  analyzed  in  Fig.  5.  Thus,  the  stratus  clouds  are  clearly 
indicated  in  this  computer  analysis  of  the  signature  slope,  even  though 
the  initial  assumption  of  a  constant  0  is  obviously  unreal.  The  presence 
of  horizon  to  1 ly-ext ensive  higher  clouds  is  suggested  by  the  areas  of 
negative  values  (signal  increase)  at  500  m.  The  hand  analysis  of  Fig. 

5  shows  these  clouds  only  in  the  high-elevation  lidar  da'a.  Thus,  the 
computer  is  capable  of  analyzing  small  changes  in  signature  slope  that 
escape  the  eye. 

A  solution  to  Eq .  (3)  for  the  optical  parameters  0  and  C  is  possible 
when  the  existing  relation  between  these  quantities  is  known.  Fenn,' 
using  the  data  of  Barteneva/'  and  Curcio  and  Kr.estrick’  have  shown  that 
a  relation  such  as 

d  In  0  ,  . .  , 

— — -  =  constant  =  k„  (o) 

d  In  r  2 


is  valid  within  20  to  30  percent  for  extinction  coefficients  between 
0.01  km  *  and  1.0  km  *  when  employing  a  broad  spectral  source.  Assuming 
such  a  relation  to  be  valid  for  the  conditions  of  this  experiment,  sub¬ 
stitution  of  Eq .  (5)  into  Eq .  (3)  yields  a  first-order  nonlinear  dif¬ 
ferential  equation: 


dR  “  C1  dR 


0 


where  =  1  4.34  k^  and  C  =  2  k^ .  The  transform  T|  =  l  g  reduces  this 
differential  equation  to  the  linear  form: 


£.  .  c  15  f 

dR  1  dR 


C 


2 


0 


(6) 
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for  which  ;  Ih-  solution  may  bo  written  as: 


(7) 


where  (_'  is  the  constant  of  integration.  The  derivative  of  S(R),  which 
is  formed  by  straight  line  segments  between  digitized  lidar- signature 
inflection  points,  is  not  unique  at  these  points — i.e.,  S(R)  is  a  piece- 
wise  di 1 1 oron t la ble  function.  Thus,  in  addition  to  the  requirement 
of  a  boundary  value,  the  integration  constant  C  must  be  re-evaluated 
for  each  S(R)  line  segment.  Solution  for  multiple  traces  requires  a 
boundary  value  for  each  trace  or  correction  and  or  assumptions  on  lidar 
variable  parameters  between  traces.  The  method  used  in  the  analysis  of 
the  present  lidar  data  was  to  obtain  a  boundary  value  of  extinction 
coefficient  from  the  slope  technique  for  part  of  a  signature  for  which 
the  backscatter  was  thought  to  bo  constant  with  range.  Solution  was 
then  taken  along  tins  trace  in  the  direction  of  the  lidar.  Boundary 
values  of  "  at  the  initial  data  point  {range  Rq)  for  each  additional 
trace  were  derived  assuming  unity  transmission  of  the  lidar  energy  to 
these  data  points  and  the  initial  value  of  C  from  the  trace  containing 
the  boundary  point.  Extinction  coefficients  were  then  evaluated  for 
each  data  point  in  the  field.  A  value  of  1.4  was  assigned  to  the  constant 
k2-  Such  a  value  probably  underestimates  the  attenuation  of  the  lidar 
energy  and  thus  increases  the  probability  of  a  valid  solution  at  regions 
of  large  signal  variations  along  the  lidar  signature.  The  validity  of 
the  solution  field  is  determined  by  re-deriving  the  S  functions  from 
the  right-hand  term  of  Eq .  (2)  from  evaluated  quantities. 


Figure  11(b)  presents  the  grid-point  analyses  of  the  a  field 
obtained.  As  in  the  previous  computer  analyses,  the  level  of  the  low 
clouds  is  clearly  shown  by  large  values  and  large  vertical  gradients 
of  - .  By  comparison  with  Fig.  10(b),  it  is  clear  that  some  of  the 
attenuation  effect  has  been  removed  from  the  data  field.  However,  the 
non-horizontal  l.soplcths  indicate  that  additional  attenuation  should  be 
accounted  for  by  employing  a  smaller  value  of  or  a  larger  value  of 
the  boundary  extinction  coefficient.  Further  experiments t ion  along 
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these  lines  was  not  considered  just  1 lied  because  ot  tile  limi.ed  quantity 
ard  quality  ol  the  data.  The  real  merit  ol  a  derived  "-field  is  that 
transmission  along  any  line  segment  doiined  by  two  points  in  the  lield 
may  be  evaluated  and  thus  slant-range  visibilities  interred. 

U.  Discussion 

Figures  10  and  11  show  the  operational  ieasibility  ol  describing 
the  spatial  distribution  of  cloud  ceiling  and  low-cloud  structure  with 
a  lidar-electronic  computer  combination.  The  results  shown  are  not 
considered  optimum  since  they  depend  on  several  assumptions  that  may 
be  relaxed  or  altered  through  further  investigations.  The  primary 
uncertainty  is  the  validity  of  the  assumed  (s-c  relation,  which  enables 
the  solution  of  Kq .  (3).  It  is  suggested  that  in  luture  experiments 
ot  the  type  described  in  this  report,  the  neodymium  laser  be  used  be¬ 
cause  of  its  broader  spectral  energy  (200  A)  in  comparison  to  that  ol 
the  ruby  laser.  The  broad  spectral  characteristic  gives  more  validity 
to  assuming  a  0-c  relation  (Twomey  and  Howell.  Ref.  5).  Also  necessary 
is  some  further  research  in  the  numerical  application  ol  Lq .  (7).  For 
example,  various  boundary  values  can  be  evaluated  by  using  a  completely 
calibrated  lidar  or  by  using  an  extinction  coefficient  derived  t tom 
the  return  of  the  lidar  energy  from  a  ground-based  target  of  known 
reflectivity  (Allen  et  al..  Ref.  6) 
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V  CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  covers  a  limited  experiment  in  the  use  of  lidar  for  the 
determination  of  ceiling  and  visibility.  The  conditions  during  the 
experiment,  although  ideal  trom  a  meteorological  point  of  view,  were  by 
no  means  optimum  in  terms  of  equipment  or  experimental  procedures.  For 
example,  the  use  of  the  neodymium  version  of  the  Mark  V  lidar,  capable 
of  iiring  every  4  to  5  seconds  would  have  been  preferable;  also  more 
complete  calibration  and  reference-measurement  procedures  are  desirable. 
However,  the  results  demonstrate  a  real  potential  of  lidar  for  use 
in  operational  determinations  of  aircraft  landing  visibility  conditions. 

It  is  evident  that  the  lidar  can  obtain  cloud  ceiling,  even  when 
the  cloud  base  is  diffuse,  at  locations  distant  from  the  point  of 
observation.  At  Hamilton  AFB  it  was  possible  to  obtain  detailed  infor¬ 
mation  on  the  cloud  conditions  at  locations  along  the  approach  path, 
where,  because  of  the  marshes  and  open  water,  conventional  ceilometers 
could  not  be  operated. 

The  possibility  of  processing  lidar  observations  to  obtain  quanti¬ 
tative  data  on  the  extinction  coefficient — i.e.,  the  optical  parameter 
significant  to  "visibility"  determinations — has  been  explored  with 
indications  that  operationally  useful  computer-produced  analyses  are 
feasible.  The  objective  determination  of  a  cross  section  showing  a 
field  of  values  that  are  related  to  atmospheric  scattering  activity 
and  that  portray  conditions  of  low-cloud  ceiling  and  reduced  horizontal 
visibility  is  considered  to  be  a  major  advance  in  this  problem  area, 
and  one  that  offers  hope  of  eventually  being  able  to  infer  "slant 
visibility"  for  aircraft  landing  operations. 

It  is  obvious,  however,  that  much  remains  to  be  done.  For  simple 
graphical  description  of  the  cloud  base,  further  work  is  necessary  on 
the  physical  interpretation  of  inflections  on  the  lidar  signature,  and 
also  on  the  technique  of  processing  and  displaying  data.  The  conversion 
of  lidar  observations  to  digital  form  for  computer  processing,  and  the 
computational  solutions  themselves  also  need  further  development. 
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Progress  in  both  these  areas  would  involve  additional  experimental 
and  development  programs.  The  programs  should  be  geared  toward  obtaining 
optimum  techniques  of  data  processing  and  display  that  can  lead  to  the 
design  of  practical  instrumentation  for  routine  operational  use. 
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